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Abstract 

The past decade has led to significant improvements in 
our understanding of the physical structure of the molec¬ 
ular cores of cold dark clouds. Observational efforts, in 
combination with improved knowledge of cloud structure, 
now provide clear evidence that the chemistry of dark 
clouds is dominated by the depletion of gaseous species 
onto grain surfaces. We outline the basis of these obser¬ 
vational efforts and show how the abundance determi¬ 
nations have moved beyond single point analyses to the 
derivation of abundance profiles. We discuss the basic 
physics of the interaction between molecules and grain 
surfaces and show that when physics is coupled into a 
chemical model there is excellent agreement, for a lim¬ 
ited set of species, between theory and observations. We 
discuss our improved understanding of cloud chemistry 
can be used as a new tool for studies of the formation of 
stars and planetary systems. 

Key words: astrochemistry - ISM: molecules - ISM: glob¬ 
ules - ISM: clouds - stars: formation 


1. Introduction 

The formation of dense molecular condensations, and even¬ 
tually stars, involves large changes in the physical proper¬ 
ties of the atomic and molecular gas. These changes also 
have very specific consequences on the chemical interac¬ 
tions of the gas and dust inside the forming cores/stars. 
In particular as the density increases molecules in the gas 
phase collide with dust grains with greater frequency and, 
if the molecules stick with any reasonable efficiency, they 
will deplete from the gas phase. Thus, the density gradi¬ 
ent that results from core condensation is accompanied 
by chemical gradients, with the inner parts of the core 
representing high density chemistry (with freeze-out onto 
grains), and the outer parts representing the original low 
density molecular composition. Indeed, these molecular 


depletions have been suggested as the primary chemical 
indicator of the earliest stages of the star formation pro¬ 
cess (Mundy & McMullin 1997; Bergin & Langer 1997). 
Such chemical structure should be particularly evident in 
dark clouds because the low temperatures (< 10 K) pre¬ 
clude thermal evaporation of the molecules frozen in the 
mantle. 

The process of star formation is by nature a dynamic 
one and molecular emission is the primary method for ob¬ 
taining information on kinematic motions inside molecu¬ 
lar clouds. Extensive efforts have been placed towards us¬ 
ing molecular tracers, typically CS and HCO + , to search 
for the presence of star-forming infall due to gravitational 
collapse (see Evans 1999; Myers, Evans, & Ohashi 2000 
and references therein). The freeze-out of molecules onto 
grain surfaces significantly reduces the effectiveness of us¬ 
ing molecular emission as a tracer of motions. Indeed the 
inside-out collapse models of Shu (1977) predict that the 
highest infall speeds are found in the exact regions where 
molecules are expected to freeze-out. 

Despite these arguments it is only recently that detec¬ 
tions of gas phase freeze-out have become commonplace. 
Until recently inherent difficulties in extracting molecu¬ 
lar abundances from gas phase emission prevented defini¬ 
tive detections of gas-phase freeze-out. Previous studies 
assumed local thermodynamic equilibrium (LTE) to esti¬ 
mate total column densities. However, molecular emission 
from tracers with high dipole moments is far from LTE. 
In addition, secondary tracers such as CO and its iso¬ 
topic variants were previously used to estimate the total 
H 2 column density, since it cannot be observed directly. 
It is now known that depletion of CO in cold regions se¬ 
riously hampers its utility as an estimator of the total 
hydrogen column density. 

The recent advent of sensitive continuum and het¬ 
erodyne arrays probing millimeter/sub-millimeter wave¬ 
lengths, along with wide-field infrared imaging devices, 
has led to an explosion of clear detections of gas phase 
depletion/freeze-out (Bacmann et al. 2002; Bergin et al. 
2002, Tafalla et al. 2002, Hotzel et al. 2002, Caselli et 
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al. 2002, Jessop & Ward-Thompson 2001, Kramer et al. 
1999, Alves et al. 1999). At this conference alone there are 
8 contributions discussing gas-phase depletions (Carey 
et al., Di Francesco et al., Feldman et al., Savva et al., 
Kontinen et al., Lai et al., Lee et al., and Peng et al., 
this volume). Observations of dust in emission or absorp¬ 
tion provide direct knowledge of the dust distribution. 
In addition, since the dust column density and mass is 
correlated with the H 2 column density and mass (Hilde¬ 
brand 1983; Gordon 1995), these observations provide the 
clearest information to date on the spatial distribution of 
H 2 molecules. Furthermore, with some geometrical as¬ 
sumptions, the radial profiles of core density can be con¬ 
structed; these methods have greatly increased our knowl¬ 
edge of the physical conditions throughout the star for¬ 
mation process (Andre et al. 2000, Alves, Lada, & Lada 
2001). Knowledge of the density and column density dis¬ 
tribution of H 2 also significantly aids the molecular ob¬ 
servations in two ways. (1) By indirectly confirming the 
location of the H 2 density and column density peak and, 
(2) by providing the density profile which helps to un¬ 
ravel the similar effects of density and abundance on ex¬ 
citation. These advances have moved chemical analyses 
beyond surveys of objects that show evidence for com¬ 
plicated chemistry (e.g. TMC-1, L134N) towards studies 
of more centrally concentrated objects that are closer to 
collapse and star formation. 

In this concise review we outline how combined stud¬ 
ies of dust emission/absorption with molecular emission 
have improved our ability to estimate molecular abun¬ 
dances and, in consequence, how our picture of dark cloud 
chemistry has been changed to one dominated by the ef¬ 
fects of freeze-out. In the following (§§2 & 3) we briefly 
outline our current theoretical understanding of gas-grain 
interactions. Furthermore we discuss key observational ef¬ 
forts that demonstrate the systematic nature of gas-phase 
freeze-out in condensed cores and how these observations 
can be qualitatively understood by a coupling of gas-grain 
chemical models to sophisticated radiative transfer mod¬ 
els. In §4 we discuss how our improved understanding 
opens the door towards using astrochemistry as a new 
and potent tool to study the process of star formation. 

2. Gas-Phase Freeze-Out: Theory 

The rate of deposition of a molecule in the gas phase onto 
a grain surface is 

k freeze—out ~ TO. vSflgr (s )■ (1) 

where a is the grain radius, v is the mean thermal velocity, 
S is the sticking coefficient, and n gr the space density of 


grains. If we use “classical” grains with a = 1000 A and 
n gr = 1.3 x 10” 12 n(IL2) cm -3 (Spitzer 1978), then the 
timescale for a molecule to freeze-out onto a grain surface 
is: 


^freeze—out 


2 x 10 9 yr 
Smn (H 2 ) 


( 2 ) 


where m is the molecular weight and S is the sticking co¬ 
efficient. CO is a typical tracer of molecular cores which 
have typical densities of n(H 2 ) ~ 10 5 cm -3 (Evans 1999). 
For a sticking coefficient of unity the freeze-out timescale 
is < 10 4 year, well above the estimated time for ther¬ 
mal evaporation (^> 10 8 yr), but below estimated ages of 
~ 1 — 10 Myr (Palla and Stahler 2000; Hartmann 2001; 
Lee and Myers 1999). Thus it is surprising that a gas- 
phase is observed at all, instead the near total freeze-out 
of the gas phase would be expected (e.g. Iglesias 1977). 
This points to the existence of some non-thermal mecha¬ 
nism to desorb or remove molecules from the grains and 
allow for the observed active gas phase chemistry. 

Various mechanisms have been identified, including 
spot heating due to cosnric-ray or X-ray impacts, chemical 
desorption, and others. The reader is referred to Williams 
(1993) for more detailed discussion. The strength of the 
molecule bond to the grain surface is a key parameter 
that determines the efficiency of evaporation. Because low 
temperatures preclude the breaking of chemical bonds, 
molecules are not expected to be chemically bound to 
grain surfaces. Instead the approaching molecule has an 
induced dipole interaction with the grain surface or man¬ 
tle and is bound through weaker Van der Waals-London 
interactions called physical adsorption (Kittel 1996). Van 
der Waals interactions are proportional to the product of 
the polarizabilities of the molecule and the nearby sur¬ 
face species (Et oc a mo iasurface)- Because of this prop¬ 
erty some species are more tightly bound to the grain 
surface than others. For example CS is expected to be 
more tightly bound to the grain surface than CO or N 2 . 
In addition, theoretical analysis has suggested that N 2 
has a weaker bond to a water ice surface than CO (Sadlej 
et al. 1995). 

Models which have achieved the best success in re¬ 
producing observations are those that incorporate these 
assumptions regarding bond strengths, resulting in a des¬ 
orption process that preferentially removes more weakly 
bound species (Bergin & Langer 1997; Charnley 1997; Ne- 
jad & Wagenblast 1999; Aikawa et al. 2001; Li et al. 2002). 
As as sample in Figure 1 we provide chemical abundances 
as a function of time for a cloud with a centrally concen¬ 
trated density profile. This result uses the model for B68 
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Log 10 H 2 Density in cm " 3 Log 10 H 9 Density in cm " 3 


4.7 5.1 5.4 5.6 5.7 4.7 5.1 5.4 5.6 5.7 



ers as a function of cloud depth (Ay = 1 to 17 mag) from 
a time dependent model sampled at t = 2 xlO 4 yr. Right 
panels: same except for t = 1 x 10 4 yr. The top axis on 
left, and right sides show the density that corresponds to 
extinction on the lower x-axis. The model shown is rep¬ 
resentative for the evolution of the chemistry in an object 
that is centrally concentrated such that the highest den¬ 
sities are present at the greatest extinctions. This gives 
rise to the observed chemical structure as the depletion 
rate increases with density. 

discussed by Bergin et al. (2002) where the density pro¬ 
file is that for a Bonnor-Ebert sphere near-equilibrium 
that is kept constant with time as the chemistry evolves. 
This is a simplification as more realistic models couple 
the physical and chemical evolution (Bergin & Langer 
1997; Aikawa et al. 2001; Li et al. 2002). Nonetheless it 
is illustrative of the observed effects. 

Figure 1 shows chemical abundances for species that 
are typically used to trace molecular clouds. At early 
times (t = 2 xlO 4 yr) the effects of freeze-out are evi¬ 
dent, but the overall effect is small. At later times (t = 7 
x 10 4 yr) strong gas-phase depletions are seen for nearly 
all species. CS has the largest abundance decrease fol¬ 
lowed by other carbon-bearing species. In contrast the 
nitrogen chemistry remains active at later times due to 
the presumed volatility of the N 2 molecule (see Bergin & 
Langer 1997). One byproduct of CO depletion is an in¬ 
crease in the DCO + /HCO + ratio (see Aikawa et al. 2001; 
Caselli et al 2002; Caselli 2002 for further discussion). In 


the following section we show that the pattern of differ¬ 
ential depletions - early CS depletion, followed by CO, 
and finally N 2 - is observed in numerous sources. Sim¬ 
ilar effects would be observed in any model that incor¬ 
porates these particular assumptions with regards to the 
gas-grain physics and evolves to a centrally concentrated 
state. 

3. Gas-Phase Freeze-Out: Observations 

Combined gas/dust investigations have resulted in the 
development of a new method to derive molecular abun¬ 
dances. This method derives abundance profiles as a func¬ 
tion of radius, as opposed to a line of sight average. More 
specifically, the observed core density profile from the 
dust observations is used as an input to one-dimensional 
(ID) radiative transfer codes (Monte-Carlo is generally 
used; e.g. Bernes 1979; Hogerheijde and van der Tak 2000). 
The temperature profile is another parameter that is con¬ 
strained by dust thermal emission or by multi-transitional 
studies of molecular emission. Variables are the abun¬ 
dance profile and velocity structure (both bulk and tur¬ 
bulent). With these constraints, and iterating over the 
variables, the radiative transfer model predicts the radial 
profile of emission which can then be compared in detail 
to the observations. This technique has been outlined by 
van Dishoeck and Hogerheijde (1999) and applications 
can be found in van der Tak et al. (2000), Tafalla et al. 
(2002), Bergin et al. (2002), and Lee et al. (2002). 

Below we discuss some recent investigations that high¬ 
light: (1) the observed depletions are differential in nature 
and are systemic in the dense ISM (Tafalla et al. 2002) 
and (2) current gas-grain chemical models are capable of 
qualitatively matching observations of a limited set of key 
species (Bergin et al. 2002). 

3.1. Tafalla et al. 2002 

In the top panels of Figure 2 we show integrated molec¬ 
ular emission maps of the centrally concentrated molec¬ 
ular core L1517B for sample species with a map of the 
1.2 mm continuum emission (Tafalla et al. 2002). The 
notable feature in this figure is the lack of correlation 
between the 1.2 mm continuum - a tracer of H 2 with 
C 18 0 and CS emission; where the continuum peaks the 
molecular emission exhibit local minima. In contrast, the 
N 2 H + emission peak is coincident with the continuum 
emission maximum. An analysis of these data using the 
physical description of the source given by the dust con¬ 
tinuum emission is shown in the lower panels (Tafalla 
et al. 2002). Here the radial profiles of emission (solid 
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L1517B 


1.2mm C 1B 0(l-0) C 17 0(1-0) N 2 H + (l-0) CS(2-1) 



100 0 -100 100 0 -100 100 0 -100 100 0 -100 100 0 -100 



Figure2. Top panels: Maps of 1.2mm continuum, C 18 0 (TO), C 17 0 (1-0), N 2 H + (1-0), and CS (2-1) emission in 
L1517B. Bottom panels: Radial profiles of continuum and integrated emissions. Observations are indicated by solid 
squares while the lines are results from models of the radiative transfer for two abundance profiles. The solid red lines 
are constant abundance models while the blue lines represent models with a rapid decrease in the central abundance. 
Taken from Tafalla et al. (2002). 


squares) are compared to model predictions. Models that 
assume constant abundance are shown as solid red lines; 
these models predict excess emission at the core center for 
the CO isotopes and CS. To reach agreement with obser¬ 
vations the abundances of these species must be sharply 
reduced in the densest regions of the core (models shown 
as solid blue lines) - presumably through freeze-out. Con¬ 
stant abundance models provide a reasonable match for 
N 2 H + , while NH 3 requires a slight abundance increase 
towards the center. 

Tafalla et al. also compare the observations directly to 
the chemical models of Bergin & Langer (1997), finding 
good agreement to models of fast-collapse. Most signifi¬ 
cantly similar results are found for all 5 cores examined in 
this work, as well as in sources observed by other studies 
(Kuiper, Langer, & Velusamy 1996; Ohashi et al. 1999; 
Lee et al. 2002; Bergin et al. 2002). This demonstrates 
the systematic nature of depletions - they are observed 
in nearly all centrally concentrated low mass cores. Thus 
differential depletions in cold dark clouds must be a chem¬ 
ical signature of the earliest stages of star foi'mation. 


3.2. Bergin et al. 2002: B68 

The dark cloud B68 is another example of a centrally 
concentrated object that exhibits emission features which 
are attributed to gas phase freeze-out. In this core the H 2 
distribution is traced through a map of dust visual ex¬ 
tinction (Ay ) derived via near-infrared extinction map¬ 
ping techniques (e.g. Alves, Lada, and Lada 2001). In 
Figure 3 we present the B68 Ay map along with a se¬ 
ries of molecular emission maps. In this core we see that 
the N 2 H + emission peaks in a shell partially surrounding 
the peak of dust extinction. Moreover, the N 2 H + peaks 
inside the much larger C 18 0 emission hole, which itself 
lies inside the CS emission depression. Analysis of these 
data differed from Tafalla et al. (2002) in that a gas-grain 
chemical model is directly linked to a radiative transfer 
model predicting the profile of integration emission. This 
allows for chemical considerations to determine how the 
abundance is structured with depth. For example, pho¬ 
todissociation limits the abundance at cloud edges and 
the competition between depletion and desorption deter¬ 
mines how the abundance declines with density. 
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FigureS. The observed emission distribution of visual extinction or total column density (color/gray scale) in the B68 
dark cloud. Contours are integrated emission maps of C 18 0 (left-hand panel), N 2 H + (middle), and CS (right-hand 
panel). These observed emission distributions show a sequence of increasing molecular depletion from CS, to C 18 0, 
and even N 2 H+ (Bergin et al. 2002; Lada et al. 2002). 


Figure 4 presents an analysis of the N 2 H + and C 18 0 
emission; in this case the open squares show the inte¬ 
grated emission as function of visual extinction. This par¬ 
ticular representation of the data is independent of cloud 
geometry, depending only on physical variables. Further¬ 
more, if the emission is in LTE and optically thin the 
integrated emission from a source with constant abun¬ 
dance would appear as a straight line in this diagram 
(e.g. Frerking et al. 1982; Lada et al. 1994). Clearly the 
C 18 0-Ay relation deviates from a straight line at high 
extinctions. Such structure in this relation could be due 
to freeze-out or high opacities. To discriminate between 
these possibilities, observations of a lesser abundant iso¬ 
tope are typically needed (and in this case confirm freeze- 
out). The dashed line is the “best fit” radial abundance 
profile predicted in the chemical model and the solid line 
is the resulting emission profile. For C 18 0 the observa¬ 
tions require significant freeze-out and also an abundance 
reduction at low extinction due to photodissociation. In 
the case of N 2 H + an abundance reduction at low extinc¬ 
tion is required, but only a small decrease in abundance 
at high Ay is needed to match observations. 

In Figure 5 we present a more qualitative comparison 
of the structure predicted by the chemical models to the 
B68 data. Here the left-hand panel represents the peak 
contours predicted by a spherical radiative transfer model 


combined with the chemistry for transitions of CS, C 18 0 
and N 2 H+. The right-hand panel shows an identical rep¬ 
resentation of the observed emission in B68. Although 
clear divergences from spherical symmetry are evident 
there exists striking agreement between model and ob¬ 
servation. The results presented here and elsewhere sug¬ 
gests that we are making progress in our understanding 
of the chemistry in these cold regions (e.g. Tafalla et al. 
2002; Aikawa et al. 2001; Li et al. 2002; Lee et al. 2002). 
The decrease in the N 2 H + abundance is also significant 
as this suggests that its parent molecule (N 2 ) is freezing 
onto grain surfaces. If this result is observed beyond B68 
it could have important implications for our ability to use 
molecules as the tracers of star formation dynamics. 

3.3. Summary of Trends 

The above results are excellent examples of the trends 
that are becoming apparent in chemical studies of low 
mass dark clouds. These trends have suggested changes in 
the way molecular emission is used to probe star forming 
cores. 

1. The use of C 18 0 as a probe of the H 2 column den¬ 
sity in regions below the sublimation temperature of 
CO (~ 20 K) is no longer appropriate. Whether this 
extends to warmer cores in giant molecular clouds re- 
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A v (mag) A v (mag) 

FigurePoint by point comparison of C 18 0 J=l-0 and N 2 H + integrated intensity with visual extinction for the 
entire B68 dark cloud. In all plots the data are presented as open squares with error bars while solid curves represent 
the emission predicted by a model combining chemistry with a Monte-Carlo radiative transfer code. The dashed lines 
are the best fit molecular abundance profiles with the axis labeled to the right (abundance of N 2 H + is multiplied 
by 100). In a symmetric model the abundance profiles represent the radial chemical structure within a pencil beam 
towards the extinction peak and thus are shown only for one-half of the total maximum extinction. Adapted from 
Bergin et al. (2002). 
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Figure 5. Comparison of model emission predictions with observational data. Left: Color scale map of CS (J=3-2) 
emission predicted from the B68 chemical model. Also shown are the peak contours of C 18 0 J=l-0 (blue) and N 2 H + 
J=l-0 (yellow). Right: Similar plot of the observed emission in B68. Although the core clearly shows deviations from 
spherical symmetry there is a close correspondence between model predictions and observations. 
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mains an open issue (see Mauersberger et al. 1992; 
Chandler & Carlstrom 1996; Gibb & Little 1998; Savva 
et al., this volume) 

2. Due to freeze-out the CS abundance can be sharply 
reduced when n > 10 4 cm -3 and is no longer a reliable 
tracer of motions in high density gas. However, this 
molecule is still a good probe of low density molecular 
gas prior to condensation and after star formation. 

3. The emission from N 2 H + and NH 3 most closely re¬ 
sembles the maps of dust emission/extinction. N 2 D + , 
DCO + , and NH 2 D are also found to be well correlated 
with the dust (Caselli et al. 2002; Shah and Wooten 
2001). These tracers are emerging as the best probes 
of the innermost regions of dense cores prior to the 
formation of a star. 

4. Once a star has been born the heating of the inner 
envelope returns some tracers to the gas phase. This 
has not specifically been addressed in this paper but 
the reader is referred to recent work in this area in¬ 
cluding: Hogerheijde et al. (1999), Maret et al. (2002), 
Jprgensen et al. (2002), and references therein. 

4. The Future of Dark Cloud Astrochemistry 

For a few key species we are now finding a convergence 
between observation and theory. This opens the door to¬ 
wards using chemistry as a new tool to study the process 
of star formation. Below we discuss how this can be ac¬ 
complished through both observational and theoretical 
efforts. 

4.1. Observations 

One new tool is the isolation of specific tracers that ei¬ 
ther do not freeze-out or freeze onto grain surfaces at late 
stages. Observations of these tracers (e.g. N 2 H + ) over 
small or large scales can provide complementary kine- 
matical information to continuum observations. For in¬ 
stance, in regions that are currently forming stellar clus¬ 
ters, we can investigated the core-core velocity dispersion 
and dense core gas motions relative to the overall struc¬ 
ture (see Walsh; Hogerheijde, this volume). On similar 
scales combined studies of an “early” depleter (e.g. CS) 
and a high density tracer (e.g. N 2 H+) can also take ad¬ 
vantage of the chemical evolution to detect cores prior 
to condensation (CS), during condensed stages (N 2 H + ) 
and after star formation (see Lee et al. 2002; Lee, this 
volume). 

The presence of differential depletions can certainly 
reduce the effectiveness of using molecular emission as a 


tracer of star formation kinematics, However, with care¬ 
ful choice of tracers one can use our expanded chemical 
knowledge to reconstruct the radial profile of kinematical 
motions. Indeed it is likely that the key to fully unravel¬ 
ing star formation dynamics lies in our improved under¬ 
standing of the chemistry. One example is in B 68 where 
the emission lines of C 18 0 and N 2 H + are characterized 
by a systematic and well defined velocity gradient across 
the face of the cloud. However, the magnitude and di¬ 
rection of the velocity gradient is different between the 
two tracers. The chemical analysis presented in Figure 4 
shows that N 2 H + emission traces the innermost regions 
of the B 68 core, while emission from C 18 0 predominantly 
traces the outer layers. Thus, if the gradient is due to ro¬ 
tation, then the combined chemical/dynamical analysis 
suggests that the core is differentially rotating (Lada et 
al. 2002 ). 

Another example of the importance of chemistry to 
dynamics is shown in Figure 8 . The left-hand panel shows 
spectra of 3 molecules taken towards the same position 
in the L183 molecular core. The double peaked structure 
observed in the optically thick CS and HCO + emission 
is now commonly used as a tracer of motions along the 
line of sight (e.g. Evans 1999). However, in L183 the CS 
profile is consistent with infalling motions, while HCO + 
is consistent with outflow. Since both of these transitions 
sample the same level of excitation this set of conflicting 
data must be due to chemistry. In particular, the answer 
must lie in the differences in the abundance structure 
of each molecule, and changes in the velocity field, as a 
function of depth. 

A qualitative example of how these tracers could sam¬ 
ple different motions along the line of sight is shown in 
the right-hand panel of Figure 8 . This example uses the 
combined chemical/radiative transfer model constructed 
for B 68 . In the plot the x-axis is the radial position inside 
the cloud with a radius of 0.065 pc, and the y-axis shows 
the time evolution. The lines delineate the r = 1 sur¬ 
face of the line core emission. All positions to the left of a 
given line for a particular transition have higher opacities. 
For example, for t < 8 x 10 4 yr emission from N 2 H 4 " is 
optically thin and traces the entire line of sight through 
the cloud. At later times the abundance rises and the 
emission becomes thick and only probes a surface that is 
successively further from the core center. 

When optically thick molecular transitions are used 
to probe motion via asymmetric self-reversals the emis¬ 
sion preferentially probes the r = 1 surface and hence the 
magnitude and direction of motions at that depth. This 
plot clearly demonstrates that as the chemistry evolves 
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Figure6. Left: Comparison of CS (2-1), HCO + (1-0), and N 2 H + (1-0) spectra taken towards the same position in the 
dark cloud L183. The vertical line shows the systemic velocity of the cloud as determined by the optically thin N 2 H + 
emission. The CS and HCO + optically thick spectra show evidence for motions along the line of sight. However, for 
CS the spectra is consistent with infall, while HCO + appears to be consistent with outflowing motions (Lee & Myers 
2002, in preparation). Right: Plot showing the emission surface for selected molecular transitions from a combined 
chemical and radiative transfer model of B68. The x-axis is the radius position in the core and the y-axis shows 
the chemical evolution. The lines delineate the r = 1 surface for the core of the spectral line of selected molecular 
transitions. All positions to the left of a given line trace greater opacities and are not probed by observations. When 
optically thick molecular line is used to trace motions via asymmetries in the profile the emission preferentially probes 
the t = 1 surface and hence the motions at that depth. This plot demonstrates that as the chemistry evolves each 
molecule is sensitive to motions at different regions along the line of sight. 


the various tracers are sensitive to different layers within 
the cloud. For instance in this model CS (J=2-l) and 
HCO + (J=l-0) initially both probe the outer layers. As 
the cloud evolves, CS differentially depletes and the r = 1 
surface retreats deeper into the cloud. Eventually, (for 
t > 3 x 10 4 yr) HCO + emission traces outer layers, while 
CS probes the inner layers. This plot is not supplied 
as a template for the actual evolution of any particu¬ 
lar core, rather it is shown as an example of how evolu¬ 
tionary changes in the chemical composition can affect 
observed emission, and more importantly the interpreta¬ 
tion. Therefore, with observation of several species, and 
knowledge of the chemical abundance and physical pro¬ 
files, another powerful new tool is the ability to search 
for structure in the radial velocity field, a key to under¬ 
standing star formation. 


5. Theory 

Our understanding of chemical theory in dark clouds also 
must evolve. Areas that clearly need further investigation 
include comparison of models to observations of similar 
cores (L1544, L1517B, B68) in species beyond the typical 
suite of tracers (e.g. CO, CS, N 2 H + ). For example Lee 
et al. (2002) find that the abundance of CCS is depleted 
in the the L1512 and L1544 cores, but find evidence for 
a small abundance rise in the center. This would be in 
agreement with the models of Ruffle et al. (1997) that sug¬ 
gest a rebound in the abundances of complex molecules 
due to CO freeze-out. Do other complex molecules pro¬ 
vide similar results and can we further refine our under¬ 
standing of gas phase pathways? This moves the testing 
of chemical models beyond comparison to abundances de¬ 
rived in template sources, such as TMC-1 and L134N, to 
observational comparisons with more centrally condensed 
objects that have well defined physical properties. Abun¬ 
dances derived in these sources still have uncertainties 
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due to calibration, collision rates, and geometrical as¬ 
sumptions, but can generally be considered more reliable 
than simple line of sight averages. Finally one challenge 
for future theory is to place the above observations in 
context of the non-detections of water vapor and molecu¬ 
lar oxygen in cold cores (Snell et al. 2000; Bergin & Snell 
2002). Some efforts have begun to examine this question 
(Viti et al 2001; Charnley et al 2001; Roberts & Herbst 
2002 ). 

Models that combine chemistry with dynamics have 
become increasingly more sophisticated. The ultimate goal 
of these efforts is to use chemical observations to refine 
our understanding of star formation theory. This goal 
is noteworthy and it remains to be seen whether the 
chemical unknowns regarding gas-grain interactions dom¬ 
inate (e.g. sticking coefficients, binding energies, grain 
sub-structure, surface chemistry, exact desorption pro¬ 
cess). Combined models have recently begun to attempt 
to reproduce not only the pattern of depletions, but also 
the column densities observed in L1544, a new template 
source (Aikawa et al. 2001; Li et al. 2002). This com¬ 
mendable step is needed to move theory beyond qual¬ 
itative towards quantitative comparisons with observa¬ 
tions. However, in §3.2 we demonstrated that we now 
have the capability to place model predictions into radia¬ 
tive transfer models and predict direct observables. This 
is important because column density estimates depend on 
assumptions regarding excitation that may not be repro¬ 
duced in the model at the given time where model col¬ 
umn densities best match observational estimates. Hence 
one could reproduce the column density and not the ob¬ 
served emission because of differences in the radial pro¬ 
file of physical conditions. This task is time consuming, 
but sophisticated radiative transfer codes are now read¬ 
ily available (see Hogerheijde & van der Tak 2002). The 
combination of chemical/dynamical models with radia¬ 
tive transfer truly offers the best opportunity to take 
chemical models and astrochemistry to the next level and 
potentially begin to test star formation theory. 

6. Summary 

The chemical processes that result in the formation of 
molecules play an important role in star and planetary 
formation through molecular contributions to the stabil¬ 
ity of the cloud against gravitational forces. Given the 
widespread nature of gas-phase freeze-out it is becoming 
increasingly evident that knowledge of the chemistry is 
also required in order to fully characterize the star for¬ 
mation process. The future of dark cloud chemical stud¬ 


ies is promising as more complex and predictive models 
will be developed that can be compared to observations 
with increasing resolution and sensitivity. This will lead 
to crucial tests of chemical theory but also presents a com¬ 
pletely new approach to star formation studies, bringing 
astrochemistry into a wider realm as a potent tool. 
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